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Abstract.  A  frequency  bandgap  is  a  range  of  wave  frequencies  that  are  prohibited  from  passing 
through  a  medium.  The  dispersion  relation,  which  links  the  frequency  to  the  wave  number,  enables 
us  to  illustrate  the  bandgaps.  In  [E.  H.  Lee,  “A  survey  of  variational  methods  for  elastic  wave 
propagation  analysis  in  composites  with  periodic  structures,”  in  Dynamics  of  Composite  Materials, 
E.  H.  Lee,  ed.,  ASME,  New  York,  1972,  pp.  122-138]  and  [E.  H.  Lee  and  W.  H.  Yang,  SIAM 
J.  Appl.  Math.,  25  (1973),  pp.  492-499]  the  dispersion  relation  was  studied  theoretically  for  the 
one-dimensional  periodic  structure  made  of  two  materials  arranged  symmetrically  with  respect  to 
the  center  of  the  cell.  Their  dispersion  relation  formulas  can  be  similarly  extended  to  a  multilayered 
symmetric  cell  configuration,  but  not  to  a  general  (nonsymmetric)  cell  configuration.  The  general 
model  was  considered  in  [M.  Shen  and  W.  Cao,  J.  Phys.  D,  33  (2000),  pp.  1150-1154],  where  each 
unit  cell  of  the  periodic  layered  structure  contains  several  sublayers  of  arbitrary  lengths  and  materials. 
Using  the  transfer  matrix  method,  the  dispersion  relation  was  successfully  derived,  involving  very 
lengthy  explicit  formulas.  In  this  paper,  we  generalize  the  work  of  Lee  and  Yang  and  develop 
recursive  dispersion  relation  formulas  for  a  general  cell  configuration.  The  recursive  formulas  are 
easy  to  implement  and,  through  several  numerical  experiments,  successfully  corroborate  the  results 
of  Shen  and  Cao. 
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1.  Introduction.  The  existence  of  bandgaps  in  one-dimensional  periodic  elas¬ 
tic  media  appears  to  have  been  first  established  by  Lord  Rayleigh  [17],  and  a  good 
review  of  the  early  work  on  wave  propagation  in  periodic  elastic  media  can  be  found 
in  Brillouin’s  classic  text  [2];  a  more  recent  comprehensive  review  that  outlines  the 
development  of  the  “band  theory”  for  electrons,  photons,  and  phonons  can  be  found 
in  Kushwaha  [10].  A  variety  of  technological  applications  has  been  suggested  for 
phononic  bandgap  materials  which  include  transducers,  acoustic  filters,  or  barriers 
for  noise  reduction,  and  even  as  a  means  for  mitigating  the  effects  of  seismic  surface 
waves.  The  study  of  phonons  and  phononic  bandgaps  associated  with  elastic  wave 
propagation  in  periodic  elastic  media  have  also  been  used  to  study  quantum  field 
effects  such  as  tunneling  phenomena  [24]. 
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Exact  dispersion  relations  for  harmonic  waves  in  an  infinite  one-dimensional 
medium  consisting  of  plane  parallel  alternating  layers  of  two  homogeneous  isotropic 
elastic  materials  were  derived  in  [20] .  Their  exact  dispersion  curves  compare  well  with 
the  “effective  stiffness”  theory  they  developed  for  the  lowest  vibrational  modes  over 
a  wide  range  of  wave  numbers.  Using  a  transfer  matrix  formalism  and  bypassing  the 
use  of  Floquet  theory,  “pseudo-”  stop  bands  and  pass  bands  were  computed  in  [3]  for 
finite,  periodically  layered  media.  They  showed  that  for  a  finite  system  containing 
at  least  ten  cells,  the  characteristics  of  the  second  stop  band  compare  well  with  that 
predicted  in  an  infinite  medium  [11],  [12]  for  a  2-2  composite  consisting  of  ceramic 
and  polymer  constituents.  This  problem  was  further  examined  in  [7],  where  it  was 
shown  that,  in  some  instances,  only  one  or  two  unit  cells  could  be  sufficient  to  depict 
the  “frequency  bandedness”  seen  in  the  infinite  medium.  Dispersion  effects  in  finite 
periodic  structures,  which  include  viscous  damping  and  the  use  of  genetic  algorithms 
for  tailoring  their  frequency  response  characteristics,  are  also  considered  in  [8]. 

Interestingly,  the  transfer  matrix  formalism  has  been  successfully  used  for  some 
time  by  geophysicists  (e.g.,  [21],  [5],  [16],  [1]),  and  for  finite,  layered  Goupillaud-type 
(equal  travel  time)  media  [4],  also  known  as  the  so-called  communication  matrix  ap¬ 
proach  [19],  [22],  [9].  These  works  are  not  usually  cited  in  prior  work  by  the  “bandgap” 
community,  but  are  included  here  to  emphasize  their  importance  in  providing  insight 
and  a  framework  for  the  analysis  of  dispersion  effects  in  periodic  elastic  media. 

Returning  our  attention  once  again  to  infinite  media,  [11]  and  [12]  study  the  dis¬ 
persion  relation  in  an  infinite  strip  of  a  periodically  repeated  cell  with  length  or  period 
a.  According  to  their  model,  the  cell  is  composed  of  two  homogeneous  elastic  mate¬ 
rials:  the  filler  (/)  and  the  matrix  material  (m).  These  materials  are  symmetrically 
arranged  with  respect  to  the  center  of  the  unit  cell  as  shown  in  Figure  1. 


m 


f 


m 


-a/2 


-b/2  0  b/2 


a/2 


x 


Fig.  1.  Symmetric  unit  cell  made  of  two  materials,  used  in  [11],  [12]. 


The  material  density  p  and  elastic  modulus  ?y  are  piecewise  constant  functions, 
taking  constant  values  with  subscripts  /  and  m  in  the  filler  and  matrix  material 
regions,  respectively.  The  density  p  and  the  elastic  modulus  r/  vary  periodically  along 
the  strip  with  position  x  and  period  a, 

f  p{x  +  a)  =  p(x), 

|  r)(x  +  a)  =  T)(x). 

Figure  2  displays  the  density  function  within  a  single  (unit)  cell,  assuming  that 
Pf  <  pm.  Since  the  cell  is  periodically  repeated  in  the  infinite  strip,  the  density  graph 
shown  in  Figure  2  is  also  periodically  repeated. 

The  general  model  was  considered  in  [18],  where  each  unit  cell  of  the  periodic  lay¬ 
ered  structure  contains  several  sublayers  of  arbitrary  lengths  and  materials.  They  were 
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-a/2  -b/2  0  b/2 


x 
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Fig.  2.  Density  in  a  symmetric  unit  cell  made  of  two  materials. 


able  to  successfully  derive  the  dispersion  relation  using  the  transfer  matrix  method, 
while  involving  very  long  explicit  formulas. 

In  this  paper,  we  generalize  the  work  of  [11],  [12]  and  develop  recursive  dispersion 
relation  formulas  for  a  general  cell  configuration.  The  recursive  formulas  are  easy  to 
implement  and,  through  several  numerical  experiments  using  Maple  [15],  successfully 
corroborate  the  results  of  [18]. 

Unlike  the  (two-material)  symmetric  cell  configuration  studied  in  [11],  [12],  the 
implicit  dispersion  relation  for  a  general  cell  configuration  appears  to  be  more  com¬ 
plex.  The  process  of  deriving  an  implicit  recursive  dispersion  relation  involves  the 
construction  of  basic  solutions  in  the  unit  cell,  which  provides  more  insight  on  how 
the  properties  of  such  solutions  relate  to  the  cell  configuration.  This  is  demonstrated 
through  two  approaches,  which  we  identify  as  the  central  expansion  approach  and  the 
quasi-symmetric  limiting  approach. 

As  shown  in  [11],  [12],  we  begin  with  the  wave  equation  with  periodic  coefficients 
rj  and  p,  described  by 


(1) 


d[v%\  _  &2U 
dx  ^  dt2 


Using  separation  of  variables,  we  assume  that  the  displacement  U(x,t)  can  be 
expressed  as 


(2) 


U(x,  t )  =  u(x)4>(t ), 


and  (1)  reduces  to  the  second  order  ordinary  differential  equation  with  periodic  coef¬ 
ficients, 


(3) 


d  du 
dx  ^  dx 


+  put2u  =  0. 


According  to  the  Floquet  theory,  for  a  fixed  ui,  the  solution  u(x )  in  (3)  is  of  the  form 
(4)  u(x)  =  v(x)eiqx, 


where  v(x)  is  a  periodic  function  with  the  same  period  a  as  the  coefficients  r)  and  p. 
Due  to  the  quasi-periodic  recursive  relation  that  follows  from  (4),  we  have 


u(x  +  a)  =  u(x)etqa, 
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Fig.  3.  Even  eigenfunction  in  a  symmetric  cell  of  two  materials  and  three  layers. 

ujx)  t 


Fig.  4.  Odd  eigenfunction  in  a  symmetric  cell  of  two  materials  and  three  layers. 


and  the  problem  of  finding  the  solution  u(x)  along  the  strip  is  reduced  to  the  single 
unit  cell  — a/2  <  x  <  a/2,  where  the  following  quasi-periodic  boundary  conditions 
apply: 


u(a/ 2)  =  u(— a/2)eiqa, 
u\a/ 2)  =  u'(—a/2)eiqa. 


The  solution  of  (3)  in  the  unit  cell,  subject  to  the  boundary  conditions  (5),  is 
then  expressed  as  a  linear  combination  of  two  eigenfunctions/linearly  independent 
solutions  of  (3).  For  convenience,  the  two  eigenfunctions  are  chosen  to  be  even  ue(x) 
and  odd  u0( x)  functions  along  the  symmetric  cell  studied  in  [11],  [12];  see  Figures  3 
and  4.  The  solution  of  (3), 


(6) 


u{x)  =  ue(x)  +  Cu0(x), 
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satisfies  the  quasi-periodic  boundary  conditions  (5)  if 


(7) 


^  _.ue{a/ 2) 

O  l  .  . 

u0(a/2) 


qa 
tan  — 
2 


and 


=  <(a/2) 

<(a/2) 


According  to  [11],  [12],  the  compatibility  of  the  two  relations  above  results  in  the 
implicit  form  of  the  dispersion  relation, 


(8) 


Ue(a/ 2)  qa  ue(a  2)  qa 

— t — —  tan  —  = - — — —  cot  — . 

uQ(a  2)  2  u'  {a  2)  2 


Further  simplifications  imply  the  following  equivalent  forms: 

,  2  (Qa\  _  u'e(a/2)u0{a/ 2) 

an  V  2  ;  u'0(a/2)ue(a/2) 

or 


(9) 


,  \  u'0(a/2)ue(a/2)  +  u'e(a/2)uo(a/2) 

qa  u'0(a/2)ue(a/2)  -  u'e{a/2)u0(a/2) 


The  dispersion  relation  is  then  obtained  after  the  construction  of  the  even  and  odd 
eigenfunctions  ue(x)  and  u0{x )  in  the  unit  cell,  and  their  substitution  into  (8)  or  (9). 
As  discussed  in  [11],  [12],  due  to  (4),  it  is  only  necessary  to  consider  the  wave  number  q 
limited  to  the  domain  0  <  q  <  n/a.  The  dispersion  relation  graph  displays  a  banded 
frequency  spectrum,  comprising  bands  which  transmit  Floquet  waves  and  no  pass 
bands  which  do  not. 


2.  Generalized  form  of  the  dispersion  relation.  In  this  section,  we  derive 
the  dispersion  relation  for  a  general  (unit)  cell  configuration,  made  of  an  arbitrary 
number  of  layers  and  materials.  The  steps  involved  are  summarized  below,  and  are 
generalizations  of  the  work  of  [11],  [12]  discussed  earlier.  Through  the  rest  of  the 
paper,  the  interval  of  the  unit  cell  with  length  a  is  [—b,d\.  Here  0  <  b,  d  <  a,  and 
b+d  =  a.  In  the  special  case  of  a  symmetric  cell  configuration  seen  before,  b=  d=  §. 

(i)  The  application  of  Floquet ’s  theorem  for  (3),  with  p  and  p  corresponding  to 
the  general  (unit)  cell  configuration,  yields  two  quasi-periodic  conditions, 

j  u(x  +  a)  =  u(x)eiqa, 

1  u'(x  +  a)  =  u\x)eiqa . 


Evaluating  the  conditions  above  at  x  =  —b,  we  obtain  the  generalization  of  the  bound¬ 
ary  conditions  (5), 


(10) 


u{d)  =  u(-b)eiqa, 
u'(d)  =  u'(-b)eiqa. 


(ii)  The  solution  u  =  u(x)  of  (3)  and  (10)  may  be  written  in  a  general  form 
similar  to  (6)  as 


(11) 


U  =  C\U\  +  C2U2- 


Here  u\  =  u\(x)  and  u2  =  u2{x)  are  two  eigenfunctions/linearly  independent  solutions 
of  (3),  to  be  constructed  later,  while  the  constants  C\  and  C2  are  unknown.  For  a 
symmetric  cell  configuration,  such  as  the  one  studied  in  [11],  [12],  u\  and  u2  become 
even  and  odd  functions. 
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(iii)  As  a  generalization  of  (7),  after  substituting  (11)  into  the  quasi-boundary 
conditions  (10),  we  obtain  a  linear  homogeneous  system  of  equations  for  the  unknowns 
Ci  and  C'2, 


(mi (d)  -  mi (-b)eiqa)  Ci  +  (u2(d)  -  u2(-b)eiqa)  C2  =  0, 
K(d)  -  u'1{—b)eiqa)  Ci  +  (u'2(d)  -  u'2(—b)eiqa)  C2  =  0. 


(iv)  Seeking  a  nontrivial/nonzero  solution,  we  set  the  determinant  of  the  sys¬ 
tem  (12)  to  zero,  a  condition  which  replaces  the  compatibility  conditions  used  earlier 
in  (7)-(8)  by  [11],  [12].  This  yields  the  w  =  ui{q)  relation  of  the  form 

(mi (dj  -  Mi(-6)e^)(M'  (d)  -  u’2{-b)eiqa) 

=  (u[(d)  -  M/1(-b)e^)(M2(d)  -  u2{—b)eiqa). 

Equation  (13)  represents  the  dispersion  relation  for  a  general  cell  configuration  in  its 
implicit  complex  form. 

(v)  After  a  few  manipulations  of  (13),  using  the  property  of  the  Wronskian 
W( x)  discussed  below,  we  obtain  the  dispersion  relation  in  its  implicit  real  form. 
The  Wronskian  W(x)  of  ui(x)  and  u2(x),  two  linearly  independent  solutions  of  the 
differential  equation  (3),  is  given  by 

W(x)  =  Ui{x)u'2{x)  —  u'1(x)u2(x). 


It  directly  follows  from  (3)  that 

-^[r}(x)W(x)]  =  0, 

and  therefore, 


W(x) 


C 

V(x)' 


where  C  is  a  constant.  Assuming  that  rj(—b)  =  rj(d),  one  deduces  that  W(—b)  =  W(d). 
Returning  our  attention  to  (13),  after  multiplying  and  reorganizing  the  terms,  we 
obtain 


(14)  W(d)  -  J(— 6,  d)eiqa  +  W(-b)ei2qa  =  0, 

where  J(—b,d)  =  ui(—b)u2(d)  +  ui(d)u2(—b)  —  Ui(— b)u2(d)  —  m,1(cZ)m2(— b). 

Under  the  assumption  that  r](—b)  =  rj(d ),  one  deduces  that  W(—b)  =  W(d)  and 
the  relation  (14)  becomes 

W(-b)[eiqa +e-iqa]  =  J(—b,  d). 

From  here  the  dispersion  relation  is  expressed  in  its  implicit  real  form  as 

J(—b,d ) 

2  =  TVFsp 

or  equivalently 


(15) 
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Further  simplifications  follow  if  u\(x)  and  U2(x)  are  chosen  to  satisfy  u\(—b)  = 
u'2(—b)  =  1  and  u[(—b)  =  U2{~b)  =  0.  Then  W(—b)  =  1  and  the  implicit  dispersion 
relation, 


2  cosfya) 


J(-6,  d) 
W(—b)  ’ 


takes  the  form 


2cos(ga)  =  u\{d)  +u'2(d). 


Similar  arguments  can  be  found  in  [14]. 

A  simplified  form  of  (15)  used  to  express  the  dispersion  relation  for  the  (two- 
material)  symmetric  cell  configuration  was  derived  in  [11],  [12]  and  given  by  (9). 
Indeed,  when  6  =  d  =  ]|,  we  obtain 


(16) 


2  cos(ga) 


m(-fK;(f)  +m(fK2(-f) 

ui(-fK(-f) 


<(-§  Wf )  -  <(§  W-f ) 

<(-§  )«2(-f) 


Due  to  the  symmetry  of  the  cell  configuration,  rti  and  112  become  even  and  odd 
functions,  respectively.  This  means  that 


f  ui(§)  =  ui(-f),  f  u2(|)  =  -u2(-f), 

l«i(§)  =  -«i(-f)  l«2(|)  =  «2(-f)- 

After  substituting  these  relations  into  (16)  and  replacing  ui  =  ue  and  112  =  ua, 
we  obtain  the  dispersion  relation  (9),  as  seen  before  in  [11],  [12]. 

As  a  conclusion,  relation  (15)  subject  to  the  boundary  requirement  r](—b)  =  rj(d) 
represents  the  generalized  form  of  the  dispersion  relation  (9)  derived  in  [11],  [12].  This 
is  also  confirmed  by  our  numerical  experiments  with  the  choice  of  the  unit  cell  boxed 
in  Figure  5,  with  border  layers  occupied  by  the  same  material,  ensuring  the  same 
value  for  ry( x)  and  therefore  the  same  value  for  W(x)  along  the  border  layers. 


Fig.  5.  Unit  cell  selection  in  a  periodic  medium  with  a  general  cell  configuration. 

3.  Recursive  formula  of  the  dispersion  relation  using  the  central  ex¬ 
pansion  approach.  In  order  to  develop  the  generalized  dispersion  relation  (15),  we 
need  to  find  two  eigenfunctions  u\{x)  and  u^{x)  of  (3)  in  a  unit  cell  of  our  choice. 
The  unit  cell  of  choice,  used  in  the  central  expansion  approach,  is  shown  in  the  last 
diagram  of  Figure  6.  This  is  obtained  after  shifting  and  renumbering  the  general  cell 
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3-Layered  Symmetric  Cell  (Lee  and  Yang,  1973) 


2  1  2 


General  Cell 


1 

2 

3 

.  .  . 

M 

Shifte 

d  Cell 

M 

1 

2 

3 

.  .  . 

M 

Renumbered  Cell 


N“ 

3- 

2- 

1 

2* 

y 

N+ 

Fig.  6.  The  stages  of  development  of  the  central  expansion  approach  for  a  given  general  cell 
configuration  of  M -layers.  Here  N+  =  N~  =  N ,  where  N  =  |~ 4f  ]  +  1  or  N  =  \ 4r~| . 


of  M-layers.  The  purpose  of  the  shifting  is  to  have  the  border  layers  made  of  the 
same  material  to  ensure  that  r](—b)  =  r/(d ),  which  implies  that  W(—b)  =  W(d)  and 
therefore  (15)  holds.  The  purpose  of  the  renumbering  is  to  create  a  central  layer  in 
the  cell  similar  to  the  three-layered  symmetric  cell  configuration  previously  studied 
in  [11],  [12].  In  the  renumbered  scheme,  N  =  |~4p]  +  1,  unless  the  cell  already  has 
an  odd  number  of  layers  M  and  the  border  layers  are  made  of  the  same  material,  in 
which  case  N  =  [4f].  The  (— )  and  (+)  superscripts  on  the  renumbered  cell  diagram 
of  Figure  6  are  used  for  the  layers  numbered  2,  3, . . . ,  IV  to  indicate,  respectively,  their 
left  and  right  positions  with  respect  to  the  central  layer.  The  central  layer  is  marked  1 
to  distinguish  it  from  the  layers  on  the  other  diagrams  marked  1. 

As  a  result,  we  may  follow  the  method  of  [11],  [12].  We  begin  with  the  solution 
in  the  central  layer  of  the  cell  and  expand  it  to  the  right  and  left  layers  using  the 
continuity  conditions  of  stress  and  displacement  at  the  layer  interfaces.  Indeed,  in 
a  given  layer  of  the  cell,  with  constant  material  properties  p  and  77  and  wave  speed 
c  =  Vd/Pi  the  solution  of  (3)  is  of  the  form 


(17) 


u(x )  =  A  cos 


+  B  sin 


In  the  notation  that  follows  p± ,  77^,  and  wave  speed  =  y  vf/pf  indicate  the  corre¬ 
sponding  values  for  the  jth  layer,  j  =  2 , . . . ,  IV,  located  to  the  right  ((+)  superscript) 
or  left  ((— )  superscript)  of  the  central  layer  with  j  =  1. 
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As  indicated  in  [11],  [12],  it  is  convenient  to  choose  the  two  eigenfunctions  u\(x) 
and  U2(x)  to  be  even  and  odd  functions  in  the  central  layer, 


(18) 

(19) 


LU 


u\(x)  =  cos  — x  )  ,  —  b\  <  x  <  d\, 

Cl 


U2(x)  =  sin  I  — x 

Cl 


LU 


—b±  <  x  <  d\. 


The  expansion  of  these  solutions  to  the  right  and  left  layers  would  be  of  the  form  (17), 
and  the  eigenfunction  ui(x)  along  the  unit  cell  [—6,  d]  can  be  given  as 

f  A^cos  -  ebv-i))  +  B f  sin  ^ r(x  -  djy-i))  for  djv-i  <  x  <  dN , 


u\{x)  =  < 


Af  cos  p(x  —  cfo))  +  B 3"  sin  ^^r( x  —  (fe)) 
Af  cos  —  di 
cos  (^x) 

A2  COS  f^r(x  +  b  1 

V  c2 


Bf  sin  ^^r(x  —  d 

B3  sin  (-^(x  +  bi 
\  c2 


A3  COS  (  -*=■  (x  +  62 
V  c3 


B3  sin  ( -^r  (x  +  62 

V  c3 


for  c?2  <  x  <  c?3, 
for  di  <  x  <  d2, 
for  —b\  <  x  <  di, 
for  — &2  <  x  <  —bi, 
for  — &3  <  x  <  — &25 


(20) 


An  cos  (jfr(x  +  bN-Sj  +  Bn  sin  ^^=r(x  +  6/v-i))  for  -&at  <x<  -&at_ 
Using  vector  notation  we  denote 


4  = 


for  j  =  1 

The  coefficients  in  (20),  derived  by  the  continuity  conditions  at  the  layer  inter¬ 
faces,  are  given  by  the  recursive  relations 


(21) 


4+i  =  m?4 > 


where  the  matrix  M~~  and  the  constant  parameters  are,  respectively,  given  by 


(22) 

and 

(23) 


cos  X~ 


zb  sin 


=F pf  sin  A f  pf  cos  A f 


,± _  1h  C.7'+l 


A+  =  %(dj  -  dj.  1),  A"  =  ^(bj  -  bj.  1),  pf  = 


W+ici 


Vi  = 


r  4  1 

r  1 1 

UJ 

— 1 

0 

,  do  —  0,  bo  —  0 
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for  j  =  1, . . . ,  N  —  1.  Here  and  are  the  coefficients  of  the  eigenfunction  u±(x) 
along  the  ?ith  right  layer,  while  A~,  B~  are  the  coefficients  along  the  nth  left  layer. 
The  layer  interfaces  are  located  at  — b  =  —Bn  <  —  6jv-i  <  •  •  •  <  —  bn  <  •  •  •  <  —b\  < 
d\  <  •  •  •  <  dn  <  ■  ■  ■  <  djv- 1  <  d/v  =  d.  Here  b\  =  d±,  b  +  d  =  a  and  bn  >  0,  dn  >  0 
for  n  =  1, 2, . . . ,  N.  Notice  that  the  eigenfunction  u\(x)  is  even  in  the  central  layer, 
but  it  does  not  necessarily  remain  even  after  it  expands  to  the  other  layers  of  the  cell; 
see  Figure  7.  The  cell  in  Figure  7  is  composed  of  five  layers.  Notice  that  at  the  layer 
interfaces,  u\(x)  develops  corners,  as  expected  from  the  stress  continuity  condition. 


Fig.  7.  Eigenfunction  u\(x)  in  a  nonsymmetric  cell  of  three  materials  and  five  layers.  Addi¬ 
tional  parameters  involved  in  (20)  are  uj  =  c\  =  =  c±  =  1,  p±  =  21  =  4,  and  p^f  =  22.  =  2. 

Similarly,  we  determine  the  eigenfunction  U2{x)  in  the  unit  cell  [—b,d\  as 

f  A] y"  cos  (^-( x  —  rfjv-i))  +  sin  (^-(a;  —  cZat_i))  for  dw-i  <  x  <  dm, 


u2(x) 


(24) 


Hg+  cos  (^r(a;  —  cZ2))  +  B^+  sin  (^r(*  —  cZ2)) 

A 2+  cos  (^r(a;  —  cZi))  +  B2+  sin  (^r(a;  ~  di)) 

sin  (^*) 

A*2~  cos  (j!={x  +  &i))  +  sin  (^r(a;  +  &i)) 

^3_  cos  +  M)  +  £3“  sin  (^(a;  +  &2)) 


cos  (^(a;  +  6m- i))  +  sin  (^-(a;  +  6jv-i)) 


for  d2  <  x  <  d3, 


for  di  <  x  <  d2, 


for  — 61  <  x  <  d\, 


for  —62  <  x  <  —61, 


for  — 6g  <  x  <  —b2, 


for  —  6/v  <  x  <  —bjsf. 


Using  vector  notation  we  denote  Vj 


for  j  =  1, 


N.  The  coefficients 
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of  (24)  are  defined  by  the  same  recursive  relations  described  in  (21), 


(25) 


uj+i  =  Mfvf, 


with 

(26) 


B*^ 


0 

1 


The  matrix  Mj  and  the  constant  parameters  do,  bo,  \J,  and  p±  for  j  =  1, . . . ,  N  —  1 
are  given  as  before  in  (22)-(23). 

Here  A*+  and  B*+  are  the  coefficients  of  the  eigenfunction  U2(x)  along  the  nth 
right  layer,  while  A*-,  are  the  coefficients  along  the  nth  left  layer.  Notice  that 
U2(x)  is  odd  in  the  central  layer,  but  it  does  not  necessarily  remain  odd  after  it 
expands  to  the  other  layers  of  the  cell;  see  Figure  8.  The  cell  in  Figure  7  is  composed 
of  five  layers.  Notice  that  at  the  layer  interfaces,  U2(x)  develops  corners,  as  expected 
from  the  continuity  of  the  stress  condition. 


Fig.  8.  Eigenfunction  U2(x)  in  a  nonsymmetric  cell  of  three  materials  and  five  layers.  Addi¬ 
tional  parameters  involved  in  (24)  are  uj  =  ci  =  cf  =  cf  =  1,  pf  =  —  =  4,  and  /)(  =  —  =  2. 


After  evaluating  u±(x)  and  U2(x)  at  the  boundaries  of  the  cell  located  at  x  =  —  b 
and  x  =  d,  and  then  substituting  these  expressions  into  the  general  implicit  formula 
of  the  dispersion  relation  (15),  we  obtain  the  following  expressions  for  the  numerator 
and  the  denominator,  respectively: 


ui(—b)u'2(d)  +  ui(d)u'2(—b)  —  U2(d)u'1(—b)  —  U2(—b)u'1(d) 

(27)  =  ^A~nB*n  +  anb*n-  -  A*+B~n  -  A*n~B+)  cos(A+  +  A^) 

+  —  (A+  A*n~  -  B*+B~n  -  A*+A*n~  +  B*n~B+ )  sin(A+  +  A^) 

CN 

and 

(28)  Ul(-b)u'2(-b)  +  U2(-b)u'1(-b)  =  —(AJ,B*n-  - 

CN 

Recall  that  our  leftmost  layer  is  made  of  the  same  material  as  our  rightmost  layer, 
and  hence  cjf  =  =  cn- 

Substituting  (27)  and  (28)  into  the  general  implicit  dispersion  relation  (15),  we 
derive  the  recursive  formula  for  the  dispersion  relation, 


RECURSIVE  DISPERSION  RELATIONS  IN  ONE-DIMENSIONAL  MEDIA 


681 


2cOs(ga)  ==  An  B^)  cos(A^  +  Aw) 

(29)  .  _A?fN~~AZB~N-  ,+  .  . 

i  (AnAn  AnAn  -\-BnBn  B  n B n  )  sin(XN-\-XN) 

A-n^n  ~An  Bn 

where  the  coefficients  are  obtained  through  the  recursive  relations  involving  (21)-(23), 
(25)-(26),  while  =  -^-(6/v  —  frjv-i)  and  A^-  =  —  d/v-i)-  In  the  numerical 

CN  CN 

experiments  included  in  section  5,  the  recursive  relation  generating  the  coefficients 
A^,B^,A*n  ,B*n  of  (29)  is  programmed  in  Maple  [15].  The  dispersion  relation 
graph  is  then  obtained  using  the  implicitplot  Maple  command  with  the  wave  number 
q  limited  to  the  domain  0  <  q  <  n /a. 

4.  Recursive  formula  of  the  dispersion  relation  using  the  quasi-sym- 
metric  limiting  approach.  This  alternative  approach  involves  the  quasi-symmetric 
cell  configuration  shown  in  Figure  9,  which  preserves  the  symmetric  arrangement  of 
the  materials  around  the  filler/central  layer,  and  allows  for  layers  of  the  same  material 
to  have  differing  lengths. 


Fig.  9.  Quasi-symmetric  configuration  of  a  unit  cell  made  of  M  materials. 


s  — >  0 


Fig.  10.  Obtaining  the  general  cell  configuration  as  a  limiting  case  of  the  quasi-symmetric  cell 
configuration. 

The  quasi-symmetric  limiting  approach  differs  from  the  central  expansion  ap¬ 
proach  in  the  way  the  dispersion  relation  is  derived  for  a  general  cell  configuration 
with  a  non-quasi-symmetric  configuration.  The  quasi-symmetric  limiting  approach 
views  the  general  cell  configuration  as  a  limiting  case  of  the  quasi-symmetric  config¬ 
uration  when  the  thickness  of  the  left  layer(s)  approaches  zero;  see  Figure  10.  In  the 
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limiting  procedure,  it  is  essential  to  start  with  a  quasi-symmetric  and  nonsymmetric 
cell  configuration.  The  same  limiting  procedure  applied  to  a  symmetric  cell  configu¬ 
ration  would  fail  to  produce  any  results,  because  in  a  symmetric  cell  configuration, 
as  the  thicknesses  of  the  left  layers  approach  zero,  so  do  the  thicknesses  of  the  cor¬ 
responding  right  layers.  Thus,  the  quasi-symmetric  nonsymmetric  cell  configuration 
is  essential,  as  it  allows  the  layer  thicknesses  to  the  left  and  right  of  the  central  layer 
to  be  arbitrary,  and  hence  independent.  This  approach  also  highlights  the  fact  that 
the  dispersion  formula  obtained  in  [11],  [12]  works  for  a  general  symmetric  cell  config¬ 
uration  but  fails  once  the  cell  configuration  becomes  quasi-symmetric.  In  summary, 
the  quasi-symmetric  configuration  is  a  critical  configuration  to  work  with,  because 
from  there  one  can  recover  the  dispersion  relation  for  the  general  cell  configuration, 
something  that  cannot  be  achieved  from  a  symmetric  configuration.  The  approach 
described  above  is  summarized  in  Figure  11  and  discussed  in  more  detail  below. 


3-Layered  Symmetric  Cell  (Lee  and  Yang,  1973) 


Fig.  11.  The  stages  of  development  of  the  quasi-symmetric  limiting  approach  for  a  general  cell 
configuration  of  M -layers. 

The  recursive  dispersion  formula  for  the  quasi-symmetric  cell  configuration  can  be 
derived  similarly  to  that  described  in  the  previous  section,  under  the  simplifications 
that  tjJ  =  rjf  =  rjj,  c~+1  =  c++1  =  cj+ 1,  pj  =  =  pj,  and  N  =  M  for  j  = 

1,2 ,...,M—  1.  With  these  simplifications,  the  recursive  dispersion  formula  (29) 
becomes 


-I ~AMB^  A+Bm  Am  B^j )  cos(Ajf+AM) 


2  cos (qa)  = 


(A~^jAm  AMAAf+B~lrBM  i! M  i! ,/  i  3 1 1 1 ('  A  ! ,  ■  A ,  , } 


(30) 
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The  coefficients  are  obtained  through  the  recursive  relations  involving  (21)-(23) 
and  (25)-(26).  Notice  that  =  f^(bM  -  bM- 1),  A^  =  (dM  -  dM- 1),  while 


A’  -  fib,  bj-i),  A+  =  f(dj  -  dj- r),  Pj  = 


for  j  =  1, . . . ,  M  —  1. 


As  for  the  general  cell  configuration,  while  the  coefficients  of  (30)  with  (+)  su¬ 
perscript  corresponding  to  the  right-hand  side  layers  remain  intact,  the  coefficients  of 
the  left-hand  side  layers  with  (— )  superscript  will  simplify  due  to  the  fact  that  the 
left-hand  side  layers  vanish;  see  Figure  10.  From  (21)-(23),  it  follows  that  the  coef¬ 
ficients  with  (— )  superscript  related  to  the  eigenfunction  u±(x)  along  the  left  layers 
are  given  by 

{A~+i  =  AJ  cos(A  J)-  Bj  sin(A~), 

BJ+ 1  =  Pj(B~  cos(A j )  +  Aj  sin(A“ )), 

Aj  =  1,  Bj  =  0,  j  =  M  —  1, . . . ,  1. 

In  the  limiting  case,  as  ( bj  —  bj- 1)  approaches  zero,  so  does  A  j  =  —  (bj  —  bj- 1)  for 
j  =  M, . . . ,  2.  As  a  result,  the  recursive  relation  (31)  becomes 

(32)  {  Am  =  Am~x  =  “  '  =  =  "  '  =  A2  > 

i  Bm  =  =  •  •  •  =  pm—iPm—2  •  PjBj  =  ■  •  ■  =  Pm—iPm—2  ■  ■  • P2B2  . 

Furthermore,  considering  that  pj  =  >  A)1"  =  A)“  =  Ai  (central  layer) ,  A,  =  1, 

and  B j  =  0,  the  relations  (32)  become 

f  -^m  =  Am~i  =  •••  =  A ■  =  ■  ■  ■  =  A2  =  cos(Ai), 


D-  _  Vm-icm  d_ 

%  —  „  „  %-i 

Pm  cm- i 


VjCM  m  cm  D-  Pi  cm  .  . 

— - S,  =  •  •  •  = - S2  = - sm(Ai). 


As  a  result,  the  coefficients  AM,  BM,  and  similarly  A)^  and  B*M  ,  simplify  to 


Am  =  cos(Ai)  and  BM  = 


Pi  cm  . 


A*m  =  -  sin(Ai)  and  B*M  = 


PiCm 


sin(Ai), 

cos(Ai). 


The  parameters  to  be  used  in  numerical  experiments  (see  Figure  10  and  Figure  14) 


do  —  0,  di—h/2,  dj  —  Zi/2  +  lh,  bo  —  0,  bj  —  1 1/2,  j  —  1, 2, ... ,  M, 


that  is, 


A+  —  Aj  —  Ai  — 


A+  =  Aj  = 


-,  A-=0,  J  =  1,2,...,  M. 


Finally,  by  substituting  (34)  and  (35)  in  the  recursive  dispersion  relation  (30)  for 
the  quasi-symmetric  configuration,  we  obtain  the  recursive  dispersion  relation  (36) 
for  the  general  configuration  of  the  unit  cell  made  of  M  layers/materials. 


SM  =  HPi^cosAl  +c°sAiB^  -A*+  ^isinAi  +  sin  AiB+ )  cos(AM) 
+  (-A^sinAi  -cosAiA^  +  bm^cosAi  -  ^  sin  Ai)  sin(AM)]- 
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Here  the  coefficients  with  a  (+)  superscript  are  generated  using  the  recursive  relations 
(21)-(23)  and  (25)-(26).  In  summary,  we  established  in  (36)  a  recursive  dispersion 
relation  for  the  general  configuration  of  the  unit  cell  made  of  M  layers/materials. 
This  formula  was  obtained  by  considering  the  quasi-symmetric  cell  configuration  as  an 
essential  intermediate  step.  For  cell  configurations  that  are  already  quasi-symmetric 
or  symmetric,  the  recursive  dispersion  relation  expressed  by  (30)  is  more  suitable,  as 
it  involves  fewer  recursive  steps. 

5.  Numerical  results:  Comparison  of  the  central  expansion  approach 
with  the  method  by  Shen  and  Cao  given  in  [18]. 

5.1.  Three- material  cell  ( M  —  3).  We  consider  a  cell  composed  of  three 
distinct  materials:  concrete,  nickel  alloy,  and  steel.  The  lengths  of  the  layers  are 
0.2  m,  0.25  m,  and  0.3  m,  respectively.  The  general  cell  diagram  in  Figure  6  illustrates 
the  unit  cell  used  in  [18]  for  M  =  3.  In  the  central  expansion  approach,  a  shifting 
and  renumbering  of  the  layers  in  the  original  cell  takes  place.  The  resulting  cell  is 
illustrated  by  the  renumbered  cell  diagram  given  in  Figure  6  with  N  =  [4^]  +  1  = 
[|]  + 1  =  3,  hence  the  need  for  a  fake  interface.  As  a  result,  the  cell  to  be  used  with  the 
central  expansion  approach  has  five  layers  and  the  materials  for  each  layer  are  steel, 
nickel,  concrete,  concrete,  and  steel.  Notice  the  introduction  of  a  fake  interface  on  the 
original  layer  of  concrete.  The  material  parameters  (elastic  modulus  and  density)  are 
given  in  the  appendix. 

The  two  graphs  of  the  dispersion  relation  obtained  using  the  central  expansion 
approach  given  in  (29),  and  Shen  and  Cao’s  formulas  in  [18],  overlap  in  Figure  12, 
demonstrating  the  consistency  between  the  two  methods.  The  recursive  relation  gen¬ 
erating  the  coefficients  ,  B*^1  of  (29)  is  programmed  in  Maple  [15].  The 

dispersion  relation  graph  is  then  obtained  using  the  implicitplot  Maple  command. 
Determining  more  accurately  the  values  of  the  circular  frequency  lo  for  a  given  value 
of  q ,  including  the  band  ends  with  qa  =  0  or  qa  =  n,  is  a  difficult  root-finding  prob¬ 
lem.  As  discussed  in  [11],  [12],  due  to  (4),  it  is  only  necessary  to  consider  the  wave 
number  q  limited  to  the  domain  0  <  q  <  it /a.  As  seen  in  Figure  12,  the  dispersion 
relation  graph  displays  a  banded  frequency  spectrum  using  the  reduced  zone  scheme 
for  the  wave  number  q ,  with  the  circular  frequency  to  in  the  ordinate.  The  banded 
frequency  spectrum  is  composed  of  pass  or  propagating  bands  and  stop  bands.  Over 
the  interval  0  <  q  <  7r/a,  bands  of  permissible  frequencies  appear,  separated  by 
forbidden  bands  (creating  bandgaps),  at  which  frequencies  no  Floquet  waves  can  be 
propagated.  This  band  structure  of  pass  and  nonpass  bands  shows  the  dispersive 
properties  of  the  medium.  Similar  comments  can  be  made  for  the  other  dispersion 
relation  graphs. 

5.2.  Four-material  cell  (M  —  4).  Similarly,  we  consider  a  cell  composed  of 
four  distinct  materials:  steel,  aluminum,  concrete,  and  nickel  alloy.  The  lengths  of  the 
layers  are  0.15  m,  0.1  m,  0.4  m,  and  0.2  m,  respectively.  As  before,  the  general  cell 
diagram  in  Figure  6  illustrates  the  unit  cell  used  in  [18],  when  M  =  4.  The  cell  used 
for  the  central  expansion  approach  is  illustrated  by  the  renumbered  cell  diagram  given 
in  Figure  6  with  N  =  [-rf]  +1=  [|]  +  1  =  3.  Notice  that  unlike  the  three-material 
case,  in  this  case  we  do  not  need  to  add  a  fake  interface  because  the  number  of  layers 
after  the  shift  is  already  odd.  The  materials  for  each  of  the  five  layers  are  nickel 
alloy,  steel,  aluminum,  concrete,  and  nickel  alloy.  We  plot  the  dispersion  relation 
using  both  methods,  as  shown  in  Figure  13.  Notice  again  how  well  the  two  graphs 
overlap. 
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Fig.  12.  Dispersion  graphs  for  a  three-material  cell:  Central  expansion  approach  versus  Shen 
and  Cao  [18],  using  the  reduced  zone  scheme  for  the  wave  number  q. 


Fig.  13.  Dispersion  graphs  for  a  four-material  cell:  Central  expansion  approach  versus  Shen 
and  Cao  [18],  using  the  reduced  zone  scheme  for  the  wave  number  q. 


6.  Numerical  results:  Quasi-symmetric  limiting  approach. 

6.1.  Comparison  of  the  quasi-symmetric  limiting  approach  with  the 
method  by  Shen  and  Cao  in  [18].  Here  we  consider  a  unit  cell  composed  of  five 
layers  arranged  in  a  general  configuration,  shown  in  Figure  14.  We  choose  aluminum 
as  material  1,  nickel  alloy  as  material  2,  concrete  as  materials  3  and  5,  and  steel  as 
material  4,  with  the  following  layer  thicknesses:  l\  =  0.1  m,  I2  =  0.05  m,  I3  =  0.4  m, 
I4  =  0.2  m,  and  I5  =  0.25  m.  The  material  parameters  (elastic  modulus  and  density) 
can  be  found  in  the  appendix.  The  two  graphs  displaying  the  dispersion  relation  using 
the  quasi-symmetric  limiting  approach  given  in  (36),  and  Shen  and  Cao’s  formulas 
in  [18],  overlap  in  Figure  15. 
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Fig.  14.  General  configuration  of  a  five-layer  cell. 


Shen  and  Cao  (2000) 
Quasi-Symmetric  Approach 


Fig.  15.  Dispersion  graphs  for  a  five- layer  cell:  Quasi-symmetric  limiting  approach  versus  Shen 
and  Cao  [18],  using  the  reduced  zone  scheme  for  the  wave  number  q. 


6.2.  Comparison  of  both  of  our  approaches  with  experimental  results: 
A  simplified  one-dimensional  model.  The  dispersion  relation  and  sound  attenu¬ 
ation  through  three-dimensional  structures  composed  of  periodically  arranged  cubic 
cells  were  experimentally  measured  in  [13].  The  cubic  single  cell  consisted  of  a  1  cm 
diameter  spherical  core  made  of  lead,  coated  with  a  2.5  mm  layer  of  silicone  rubber. 
The  coated  spheres  were  periodically  arranged  in  a  8  x  8  x  8  cubic  crystal  with  lattice 
constant  of  1.55  cm,  and  with  epoxy  as  the  surrounding  matrix  material.  The  cross 
section  of  the  cell  is  displayed  on  the  upper  part  of  Figure  16.  The  one-dimensional, 
three-material  symmetric  cell  model,  shown  in  the  lower  part  of  Figure  16,  may  be 
viewed,  as  suggested  by  Wang  et  al.  [23],  as  a  simplified  one-dimensional  counterpart 
of  the  three-dimensional  structure  studied  in  [13].  In  our  one-dimensional  model, 
lead  is  considered  as  material  1,  silicone  rubber  as  material  2,  and  epoxy  as  mate¬ 
rial  3  with  the  following  layer  lengths:  1  cm  (central),  with  0.25  cm  and  0.025  cm 
on  each  side.  The  material  parameters  (elastic  modulus  and  density)  are  included  in 
the  appendix.  Our  intent  here  is  to  illustrate  the  fact  that  our  dispersion  relations 
correctly  predict  bandgaps  due  to  a  so-called  localized  resonance  phenomenon  that 
has  been  observed  in  three-dimensional  ternary  systems  [13];  the  localized  resonance 
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Fig.  16.  One- dimensional  unit  cell  motivated  from  the  three-dimensional  cubic  structure  (cross 
section)  studied  in  [13]. 


can  reduce  the  magnitude  of  the  bandgap  by  two  orders  of  magnitude  relative  to  that 
caused  by  Bragg  scattering,  and  this  observation  has  spurred  renewed  interest  by  the 
acoustic  bandgap  community  for  design  of  acoustic  attenuators.  In  [6],  Bragg’s  law 
f  =  n*  v/2  a,  for  example,  predicts  a  scattering  frequency  of  about  62.5  KHz  using 
steel  scatterers  embedded  in  an  epoxy  matrix  with  a  lattice  constant  of  20  mm.  Here, 
v  is  the  longitudinal  wave  speed  of  the  matrix,  and  a  is  the  distance  between  the 
centers  of  the  scatterers  (n  =  1,  2, 3, . . .). 

Liu  et  al.  [13]  measured  acoustic  transmission  T  as  a  function  of  frequency  (250  Hz 
to  1600  Hz)  by  placing  a  receiving  transducer  at  the  center  of  their  sonic  crystals  with 
an  external  sound  source.  The  lowest  values  of  T  correspond  to  wave  frequencies 
that  are  attenuated  by  the  structure,  whereas  the  highest  values  of  T  correspond  to 
wave  frequencies  that  easily  propagate  throughout  the  structure.  Their  experiments 
reveal  that  peak  transmission  frequencies  are  located  at  /  =  600  Hz  and  /  =  1600  Hz. 
Between  these  two  frequencies,  the  transmission  coefficient  is  low.  The  few  measure¬ 
ments  made  for  low  frequency  (/  ~  300  Hz)  suggest  that  waves  with  low  frequencies 
easily  propagate  through  the  structure. 

Due  to  the  symmetric  cell  configuration  of  the  one-dimensional  simplified  model 
given  in  Figure  16,  the  central  expansion  approach  and  the  quasi-symmetric  limiting 
approach  share  the  same  dispersion  relation  formula  (29).  The  graph  in  Figure  17 
displays  the  dispersion  relation  predicted  by  (29)  and  obtained  using  the  implicitplot 
Maple  command.  Here  we  use  frequency  /  on  the  ordinate  instead  of  the  circular 
frequency  uj.  Determining  more  accurately  the  values  of  the  frequency  /  for  a  given 
value  of  q ,  including  the  band  ends  with  qa  =  0  or  qa  =  n,  is  a  difficult  root-finding 
problem.  As  seen  in  Figure  17,  between  the  frequency  range  of  0-2000  Hz,  the  graph 
exhibits  a  pass  band  for  low  frequencies  under  140  Hz,  and  what  appear  to  be  four 
additional  narrow  pass  bands  centered  approximately  at  /  =  873,  /  =  950,  /  =  1353, 
and  /  =  1907  Hz.  Closer  inspection  of  an  expanded  view  of  the  second  pass  band 
shows  that  it  is  centered  at  f  —  873.5  Hz  (Figure  18).  Expanded  views  of  the  third 
and  fourth  bands  (not  shown  here)  located  at  approximately  950  Hz  and  1353  Hz 
are  essentially  flat  to  within  numerical  roundoff,  hence  the  group  velocity  at  these 
frequencies  is  zero,  i.e. ,  vg  =  ^  =0. 

In  conclusion,  our  one-dimensional  model  appears  to  qualitatively  predict  the 
acoustic  response  of  the  three-dimensional  ternary  structure,  with  the  two  narrow  pass 
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Fig.  17.  Graph  of  the  dispersion  relation  for  one- dimensional  symmetrical  cell  model  motivated 
from  the  three-dimensional  model  of  [13],  using  the  reduced  zone  scheme  for  the  wave  number  q. 


Fig.  18.  Expanded  view  of  the  second  pass  band  illustrated  in  Figure  17. 


bands  corresponding  to  the  peaks  and  the  two  large  bandgaps  corresponding  to  the 
ranges  of  frequencies  that  are  highly  attenuated.  Our  numerical  simulations  predict 
bandgaps  with  a  lattice  constant  two  orders  of  magnitude  smaller  than  the  relevant 
wavelength,  suggesting  a  so-called  localized  resonance  phenomenon  already  observed 
in  three-dimensional  ternary  systems  [13].  However,  the  locations  of  the  bandgaps  do 
not  match  those  that  are  experimentally  observed.  This  is  to  be  expected,  considering 
the  fact  that  our  band  structure  equations  are  derived  for  infinite  one-dimensional 
periodic  elastic  media  whereas  the  experiments  in  [13]  were  conducted  on  a  finite 
three-dimensional  structure,  with  bandgaps  only  partially  developed. 

Appendix.  The  following  are  the  elastic  modulus  p  and  density  p  of  the  materials 
selected  for  the  numerical  experiments: 

1.  Concrete:  p  =  33  •  109  Pa  and  p  =  2400  kg/m3. 

2.  Steel:  p  =  210  •  109  Pa  and  p  =  7800  kg/m3. 

3.  Aluminum:  77  =  69  -  109  Pa  and  p  =  2710  kg/m3. 
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4.  Nickel  Alloy:  77  =  214  •  109  Pa  and  p  =  8130  kg/m3. 

5.  Lead:  77  =  40.8  •  109  Pa  and  p  =  11600  kg/m3. 

6.  Silicone  rubber:  77  =  117500  Pa  and  p  =  1300  kg/m3. 

7.  Epoxy:  77  =  4.4  •  109  Pa  and  p  =  1180  kg/m3. 

Acknowledgment.  The  authors  would  like  to  acknowledge  the  two  reviewers  of 
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A  ZUREK 

F  ADDESSIO 

PO  BOX  1663 

LOS  ALAMOS  NM  87545 

7  DIRECTOR 

SANDIA  NATL  LABS 
J  BISHOP  MS  0346 
E  SHERTEL  JR  MS  0382 
W  REINHART  MS  1181 
T  VOGLERMS  1181 
L  CHHABILDAS  MS  181 1 
M  FURNISH  MS  1168 
M  KIPP  MS  0378 
PO  BOX  5800 

ALBUQUERQUE  NM  87185-0307 

1  DIRECTOR 
LLNL 

M  J  MURPHY 
PO  BOX  808 
LIVERMORE  CA  94550 

3  CALTECH 

M  ORTIZ  MS  105  50 
G  RAVICHANDRAN 
T  J  AHRENS  MS  252  21 
1201  E  CALIFORNIA  BLVD 
PASADENA  CA  91125 

5  SOUTHWEST  RSRCH  INST 
C  ANDERSON 
K  DANNEMANN 
T  HOLMQUIST 
G  JOHNSON 
J  WALKER 
PO  DRAWER  28510 
SAN  ANTONIO  TX  78284 

1  TEXAS  A&M  UNIV 

DEPT  OF  MATHEMATICS 
J  WALTON 

COLLEGE  STATION  TX  77843 


2  SRI  INTERNATIONAL 
D CURRAN 

D  SHOCKEY 

333  RAVENS  WOOD  AVE 

MENLO  PARK  CA  94025 

1  VIRGINIA  POLYTECHNIC  INST 
COLLEGE  OF  ENGRG 
RBATRA 

BLACKSBURG  VA  24061-0219 

8  UNIV  OF  NEBRASKA 
DEPT  OF  ENGRG  MECH 
D  ALLEN 
F  BOBARU 
Y  DZENIS 
G  GOGOS 
M  NEGAHBAN 
R  FENG 
J  TURNER 
Z  ZHANG 
LINCOLN  NE  68588 

1  JOHNS  HOPKINS  UNIV 
DEPT  OF  MECH  ENGRG 
K  T  RAMESH 
LATROBE  122 
BALTIMORE  MD  21218 

1  WORCESTER  POLYTECHNIC  INST 
MATHEMATICAL  SCI 
K  LURIE 

WORCESTER  MA  01609 

4  UNIV  OF  UTAH 
DEPT  OF  MATH 
A  CHERKAEV 
E  CHERKAEV 
E  S  FOLIAS 
R  BRANNON 

SALT  LAKE  CITY  UT  841 12 

1  PENN  STATE  UNIV 

DEPT  OF  ENGRG  SCI  &  MECH 
F  COSTANZO 

UNIVERSITY  PARK  PA  168023 

3  UNIV  OF  DELAWARE 
DEPT  OF  MECH  ENGRG 
T  BUCHANAN 

T  W  CHOU 
M  SANTARE 
126  SPENCER  LAB 
NEWARK  DE  19716 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  UNIV  OF  DELAWARE 

CTR  FOR  COMP  ST  MATRLS 
J  GILLESPIE 
NEWARK  DE  19716 

1  COMPUTATIONAL  MECH 
CONSULTANTS 
J  A  ZUKAS 
POBOX  11314 

BALTIMORE  MD  21239-0314 

1  LOUISIANA  STATE  UNIV 
R  LIPTON 

304  LOCKETT  HALL 
BATON  ROUGE  LA  70803-4918 

2  INST  OF  ADVANCED  TECH 
UNIV  OF  TX  AUSTIN 

S  BLESS 
H  FAIR 

3925  W  BRAKER  LN  STE  400 
AUSTIN  TX  78759-5316 

1  APPLIED  RSCH  ASSOCIATES 
D  E  GRADY 

4300  SAN  MATEO  BLVD  NE 
STE  A220 

ALBUQUERQUE  NM  871 10 

1  INTERNATIONAL  RSRCH 
ASSOC  INC 
D  L  ORPHAL 
4450  BLACK  AVE 
PLEASANTON  CA  94566 

1  AKT  MISSION  RSRCH  CORP 
M  EL  RAHEB 
23052  ALCALDE  DR 
LAGUNA  HILLS  CA  92653 

1  UNIV  OF  ILLINOIS 

DEPT  OF  MECHL  SCI  &  ENGRG 
A  F  VAKAKIS 
1206  W  GREEN  ST  MC  244 
URBANA  CHAMPAIGN  IL  61801 

1  UNIV  OF  ILLINOIS 
ARSPC  ENGRG 
J  LAMBROS 

104  S  WRIGHT  ST  MC  236 
URBANA  CHAMPAIGN  IL  61801 


2  WASHINGTON  ST  UNIV 
INST  OF  SHOCK  PHYSICS 
Y  M  GUPTA 
J  ASAY 

PULLMAN  WA  99164-2814 

1  ARIZONA  STATE  UNIV 

MECHCL  &  ARSPC  ENGRG 
D  KRAJCINOVIC 
TEMPE  AZ  85287-6106 

1  NORTHWESTERN  UNIV 

DEPT  OF  CIVIL  &  ENVIRON  ENGRG 
Z  BAZANT 

2145  SHERIDAN  RD  A135 
EVANSTON  IL  60208-3109 

1  UNIV  OF  DAYTON 
RSRCH  INST 

N  SBRAR 

300  COLLEGE  PARK 
MS  SPC  1911 
DAYTON  OH  45469 

2  TEXAS  A&M  UNIV 

DEPT  OF  GEOPHYSICS  MS  3115 
F  CHESTER 
T  GANGI 

COLLEGE  STATION  TX  778431 

1  UNIV  OF  SAN  DIEGO 

DEPT  OF  MATH  &  CMPTR  SCI 
AVELO 

5998  ALCALA  PARK 
SAN  DIEGO  CA  92110 

1  NATIONAL  INST  OF 

STANDARDS  &  TECHLGY 
BLDG  &  FIRE  RSRCH  LAB 
J  MAIN 

100  BUREAU  DR  MS  8611 
GAITHERSBURG  MD  20899-861 1 

1  MIT 

DEPT  ARNTCS  ASTRNTCS 
R  RADOVITZKY 
77  MASSACHUSETTS  AVE 
CAMBRIDGE  MA  02139 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  MIT 

DEPT  MATLS  SCI  ENGRG 
E  THOMAS 

77  MASSACHUSETTS  AVE 
CAMBRIDGE  MA  02139 

2  MATERIALS  SCI  CORP 
A  CAIAZZO 

R  LAVERTY 

181  GIBRALTAR  RD 

HORSHAM  PA  19044 

2  DIR  USARL 
AMSRD  ARL  D 
C  CHABALOWSKI 
V  WEISS 

2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

ABERDEEN  PROVING  GROUND 

81  DIR  USARL 

AMSRD  ARL  WM 
SKARNA 
J  MCCAULEY 
P  PLOSTINS 
J  SMITH 
T  WRIGHT 
AMSRD  ARL  WM  B 
J  NEWILL 
M  ZOLTOSKI 
AMSRD  ARL  WM  BA 
D  LYON 

AMSRD  ARL  WM  BC 
P  WE1NACHT 
AMSRD  ARL  WM  BD 
P  CONROY 
B  FORCH 

R  PESCE  RODRIGUEZ 
BRICE 

AMSRD  ARL  WM  BF 
W  OBERLE 
AMSRD  ARL  WM  M 
R  DOWDING 
S  MCKNIGHT 
AMSRD  ARL  WM  MA 
J  ANDZELM 
R  JENSEN 
A  RAWLETT 
M  VANLANDINGHAM 
E  WETZEL 
AMSRD  ARL  WM  MB 
M  BERMAN 
T  BOGETTI 


M  CHOWDHURY 
W  DE  ROSSET 
W  DRYSDALE 
A  FRYDMAN 
D  HOPKINS 
L  KECSKES 
T  H  LI 

S  MATHAUDHU 
M  MINNICINO 
B  POWERS 
J  TZENG 

AMSRD  ARL  WM  MC 
R  BOSSOL1 
S  CORNELISON 
M  MAHER 
W  SPURGEON 
AMSRD  ARL  WM  MD 
J  ADAMS 
B  CHEESEMAN 
ECHIN 
KCHO 
B  DOOLEY 
C  FOUNTZOULAS 
G  GAZONAS 
J  LASALVIA 
P  PATEL 
C  RANDOW 
J  SANDS 
B  SCOTT 
CF  YEN 

AMSRD  ARL  WM  SG 
T  ROSENBERGER 
AMSRD  ARL  WM  T 
P  BAKER 
M  BURKINS 
AMSRD  ARL  WM  TA 
S  SCHOENFELD 
AMSRD  ARL  WM  TB 
N  ELDREDGE 
J  STARKENBERG 
AMSRD  ARL  WM  TC 
T  BJERKE 
T  FARRAND 
K  KIMSEY 
M  FERMEN  COKER 
D  SCHEFFLER 
S  SCHRAML 
S  SEGLETES 
AMSRD  ARL  WM  TD 
S  BILYK 
D  CASEM 
J  CLAYTON 
D  DANDEKAR 
N  GNIAZDOWSKI 


NO.  OF 

COPIES  ORGANIZATION 


M  GREENFIELD 
R  KRAFT 
B  LOVE 

M  RAFTENBERG 
E  RAPACKI 
M  SCHEIDLER 
T  WEERASOORIY A 
AMSRD  ARL  WM  TE 
J  POWELL 
B  RINGERS 
G  THOMSON 
AMSRD  ARL  WM  UV 
S  WILKERSON 
AMSRD  ARL  VT  RP 
J  BORNSTEIN 


